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GERADOR DE IMPULSOS COM OCORRÊNCIA DE POISSON 


RESUMO 


Este gerador baseia-se numa fonte radio-activa e uti- 
lisa um detector de elevada eficiência no sentido de se 
obter uma distribuição de Poisson. A verificação do 
funcionamento estatístico do gerador é indicada de forma 
sucinta. Para além da sua simplicidade é de notar que 
o gerador se constrói fácilmente com equipamento labo- 
ratorial usual, 


Os geradores de sinais com características 
aleatórias têm uma utilização cada vez mais vasta 
nos domínios da tecnologia e da investigação. 
Entre estes geradores iremos considerar os que 
produzem sequências de impulsos com ocorrência 
aleatória; se o número de impulsos por inter- 
valo de tempo, seguir uma distribuição de Pois- 
son, a sequência tem propriedades importantes 
para o processamento e análise de sinais. 

Compreende-se, assim, a importância de saber 
como obter um gerador de impulsos cuja ocor- 
rência tenha uma estatística que se aproxime da 
lei de Poisson. Esta nota descreve um esquema 
muito simples e resume os ensaios efectuados 
para verificação das propriedades do gerador. 
Mostra-se que o esquema se socorre de apare- 
lhagem normalmente existente em todos os labo- 
ratórios e que o gerador cobre a banda das audio- 
-frequências com uma estatística de Poisson cujo 
elevado nível de confiança se menciona. 


Eng.º GOMES CERVEIRA 
Assistente da C.E.E.N. e do L.S.T. (*) 


Eng.º ONOFRE MOREIRA 
Assistente da C,E,E.N, 


SYNOPSIS 


It refers to a random pulse generator based on a 
radioactive source where the detector has high efficiency 
in order to achieve a Poisson distribution as close as 
possible. Tests on the pulse distribution are briefly 
described. Besides its simplicity it can be easily assembled 
with electronic standard equipment, 


1— DESCRIÇÃO DO GERADOR 


O esquema utilizado baseia-se numa fonte de 
partículas x cuja distribuição se sabe ser bas- 
tante próxima da de Poisson. Para não afectar a 
distribuição de ocorrência das partículas « deve 
utilizar-se um detector de elevada eficiência, isto 
é, que reaja, qualquer que seja a energia das par- 
tículas «, e escolheu-se, portanto, um detector do 
tipo semicondutor. 

Como emissor de partículas « utilizámos uma 
fonte de referência AMR3 com 3.10” des./min. 
e como detector utilizámos o TME] 150-100 
fabricado pela Ortec. Quer a fonte radio-activa, 
quer o detector semicondutor se colocaram dentro 
de uma câmara de latão, a cerca de 17 mm um 
do outro. 

Na fig. 1 indica-se sucintamente a constituição 
da câmara, cujo fim é o de evitar os perigos de 
contaminação radioactiva por parte do opera- 


(*) Em comissão de serviço na Universidade de Luanda. 
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s s 
U)— FONTE AMR3 COM Ixl0 des/... 


(2)— DETECTOR 


FMEJ 150 — 100 


Fig. 1 


dor. Por meio de um contador Geiger verificou-se 
que este objectivo se encontrava perfeitamente 
atingido. Além disso a câmara foi projectada de 
forma a que fosse possível evacuá-la e traba- 
lhar-se a pressões inferiores à atmosférica. Todos 
os resultados aqui apresentados foram, todavia, 
obtidos à pressão atmosférica normal, 

A aparelhagem electrónica necessária inclui 
uma fonte de alimentação para excitação do 
detector. Na cadeia de processamento é necessá- 
rio recorrer, tal como se indica na fig. 2, a um 
pre-amplificador, a um amplificador linear e a um 
gerador de impulsos comandado (univibrador). 


Fig. 2 


Esta última unidade permite obter impulsos com 
uma forma padronizada e independente da forma 


222 


e amplitude do impulso que lhe da origem. 
Incluiu-se igualmente uma escala binária, o que 
permitiu obter uma forma de onda telegráfica 
aleatória cuja finalidade foi a de verificar o fun- 
cionamento estatístico do gerador. Este equipa- 
mento foi totalmente construído no Centro de 
Estudos de Electrónica mas, dada a gama de 
frequências de trabalho, as suas características 
não são de modo nenhum especiais ou cri- 
ticas. 

As formas de onda correspondentes à sequên- 
cia de impulsos e à sequência telegráfica, tal como 
se obtiveram num osciloscópio Tektronix 545 são 
indicadas na fig. 3. O valor de pico dos impulsos 
é de 10 Volt e a frequência média de ocorrência 


4 


é de cerca de 125 impulsos por segundo. 


* 


a) 


' + + ” + . 
, 


a) Onda de impulsos: vertical 2v/cm, horizontal 
zo ms/cm (a linha de base foi suprimida por 


razões fotográficas; 


b) Onda telegráfica: vertical 2v/cm, horizontal 


20 ms/cm., 


Fig. 3 
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[I— VERIFICAÇÃO DO FUNCIONAMENTO 
ESTATÍSTICO DO GERADOR 


Os testes para verificação do funcionamento 
estatístico do gerador foram de três tipos: 


a) Sobre a distribuição estatística do número 
de impulsos por unidade de tempo; 

b) Sobre a distribuição estatística da duração 
dos intervalos de tempo entre impulsos; 

c) Sobre a distribuição espectral da onda 
telegráfica. 


Os testes referidos nas alíneas a) e b) baseiam-se 
sobre a sequência de impulsos, mas no teste a 
que se refere a alínea c) foi utilizada a sequência 
telegráfica. 

Inicialmente utilizaram-se nas alíneas a) e b) 
testes de hipóteses baseados na máxima verosi- 
milhança [1] e no a de Pearson [1]. À exposição 
teórica e a aplicação experimental ao gerador 
destes testes já se encontra publicado [4] bem 
como a conclusão de aceitar com elevado grau 
de confiança a hipótese de o gerador ter estatis- 
tica de Poisson. 

Posteriormente, recorreu-se aos testes de hipó- 
teses baseados nas estatísticas de Cramer-von 
Mises [2] e de Kolmogorov-Smirnov [2] cujos 
resultados vieram confirmar a aceitação da 
estatística de Poisson. 

A verificação mencionada na alínea c) baseia-se 
em que, se numa onda telegráfica de valores + E 
as mudanças de estado ocorrerem segundo uma 
distribuição de Poisson de parâmetro v, o Espectro 
Quadrático W (f) tem a seguinte expressão [3], 
excepto para a frequência nula: 


W (£) = [RMS (£)]º = 

E amo 

7. (2u)) + (2x8) 

em que por RMS (f) se representa o espectro 
do valor eficaz da tensão. 


Considerando valores de frequência f >> u/= = 
40,2 Hz verifica-se que 


Log RMS (f) = - Log a 


— Log | (24 - (2=8)'] -K, — Log f 
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isto é, Log RMS (f) é aproximadamente linear 
em logaritmo de f. Como da expressão anterior 
se tira que 


RMS (f) f 


Bo RMS (E) CBF, 


e porque, por definição 


»0 Log RMS (B | RMS (f) 
“81 RMS (F,) RMS (£,) lap 


f É A 
Log, — = | >— 
52 Fo Elo 


conclui-se, finalmente, que se se medir RMS (f) 
em decibel e f em oitava, então a inclinação da 
recta é aproximadamente — 6 db/oit. 


“RMS (f) E RMS(£ 
| RMS (f,) | ada? 2 (E) 
— 20. Log = — 20.Logy 2.Log, é 


Utilizou-se este facto para verificar o funcio- 
namento do gerador, pois alterações em relação 
a ele indicar-nos-iam afastamentos em relação 
à lei de Poisson sobre a qual se baseia. Deter- 
minou-se para tanto, a densidade espectral 
do valor eficaz RMS (f) para vários valores 
de frequência e comparam-se os resultados expe- 
rimentais com aquele facto teórico. 


4 
; | -10 Ee a 
EMso -15 
RMS(D| db . Nas 
-25 5 
-30 à 
"35 Inclinação=5.967 db, 


a 
no 
tava E 


O4 08 16 315 63 125 25 5 10 20 


f (kHz) —e 


Fig. 4 
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Na determinação experimental da densidade 
espectral do valor eficaz, recorreu-se a um ana- 
lisador espectral (2107 Briúel & Kjaer Spectro- 
meter) e o registo dos valores assim obtidos foi 
feito num registador de valor quadrático médio 
(2305 Brúel & Kjaer Level Recorder). Construído 
o gráfico do nível quadrático, em db utilizando 
uma escala logarítmica para a frequência, obte- 
ve-se o resultado patente na fig. 4. 

A linha recta que melhor se adapta aos valores 
experimentais (a cada um dos quais se associou 
uma indeterminação de 1 db) tem uma inclina- 
ção de-5.967 dby/oit., valor que concorda favo- 
râvelmente com - 6 db oit., valor teórico da incli- 
nação. 

O afastamento em relação a esta recta só é 
sensível fora dó intervalo 


(100 Hz, 16 kHz) 


pelo que neste intervalo que se verifica ser a 
banda das audio-frequências, o gerador de impul- 
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sos se pode considerar como tendo ocorrência de 
Poisson. Este teste permitiu, pois, definir uma 
banda de frequências de utilização do gerador 
para a qual se pode garantir que a estatística se 
aproxima da distribuição de Poisson. 
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DRENAGEM DE 


RESUMO 


4 solução do problema de drenagem de um fluido 
Newtoniano numa placa infinita é bem conhecida [1). 
U mesmo problema aplicado a fluidos que seguem 08 mo- 
delos biparameétricos de Bingham e de Ostwald-de Waele 
não é do conhecimento do autor. Neste artigo pretende-se 
preencher tal lacuna. 


1. INTRODUÇÃO 


O problema do escoamento de uma camada 
residual de fluido newtoniano incompressível 
numa placa infinita é tratado [1] mediante a sub- 
stituição da equação da velocidade média para 
o escoamento em película descendente, em regime 
laminar estacionário, com linhas de corrente 
rectilíneas, na equação de balanço material rea- 
lizado sobre o sistema, 

O trabalho presente pretende estender o resul- 
tado a fluidos incompressíveis que seguem mo- 
delos biparamétricos (modelo de Bingham e de 
Ostwald-de Waele), com a utilização das mesmas 
hipóteses, referentes à velocidade. 


2. MODELO DE BINGHAM 
Este modelo segue a lei: 


dv, 


az — — Éo dx * Co, para | Cuz | 2% (1) 


= para |Z,,| TT, (2) 


Dado o sistema de eixos coordenados que será 
utilizado, Z,, é negativo e portanto a equação (1) 
transforma-se em: 


— — T (1º) 


—— — —  — —— mo na ———— — 0 e ea 


C. D. U. 621.603 


FLUIDOS NÃO-NEWTONIANOS 


por M. MARGARIDO RIBEIRO (*) 


SYNOPSIS 


The problem of a Newtonian fluid drainage on an 
infinite flat surface has a well known solution [1). The 
solution for the same problem when applied to fuids 
wich follow hyperparametrical modeis (Bingham's and 
Ostwald-de Waele's models) is unknown by the author. 
This paber is a trial to solve that problem, 


Para o sistema considerado (fig. 1), a lei de dis- 
tribuição das tensões de corte é: 


X 
NIVEL INICIAL DE LIQUIDO 


NIVEL AO FIM DO 
TEMPO 


Fig. 1 — Perfil de escoamento, numa secção 
normal à placa. 


Tu =tg cos É (x —3) (3) 


e portanto, o escoamento tampão começa a ma- 
nifestar-se nos pontos : 


0% sum: (3 cnc (4) 


(*) Assistente do Curso de Engenharia Química da Universidade de Luanda. 
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Daqui resulta que a espessura da camada em 


“ o 
escoamento tampão é constante e igual a A + 
é & COS 


2.1. Balanço material 


A realização de um balanço material sobre o 
elemento de volume tracejado (fig. 1), de com- 
primento unitário segundo o eixo dos yy con- 
duz a: 


-|| <v> —: | “z A f- 
Z+ AZ à dão aa: ai 


Dividindo (5) por dz Àt e passando ao limite 
vem: 


. 


d l dv do a o 9% 

(Cv, >) + E (5) =— 
dz ] di" dz dt 
(6) 


( dd 00" 
notar que — = 
dt dt 


2.2. Determinação das expressões das velo- 
cidades 


2.2.1. Velocidade média na zona de gradiente não nulo 


A substituição de (1') em (3) seguida de expli- 


dv 
* conduz a: 


citação em ordem a ” 
x 


dv, 


dx 


Ro o 


Integrando e obrigando à condição fronteira 
de ser v,==0 para x=0: 


+ 


ç 2g cos É Em =) (AR " 
— —— | 0% e —— | não” atm 9% 8 
z 
Po 2 / Ro 
e portanto, como: 
[ 
v, dx 


* o PÇ o 
Co Vga e 


(9) 
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es. 
vira: 


«2 
— pgcosb 5º 


10 
EV Ra 3 (10) 
2.22. Velocidade na zona de gradiente nulo 
og cos É 5% 
v=|v,| = ea (11) 
s Po 2 
2.3. Determinação do valor de à 
E xá sg cos E 3% 
<v, > 5! = tE TO (12) 
eo 3 
E . 
d a pg cos É 9º do 
imo (E DA) SS — 13 
| z ( “e Lo dz ( ) 
d cos b 5º 
pie a LES (14) 
d q” Es 


Substituindo as equações (4), (13) e (14) em 
(6), vem: 


[E 24 cos (8 92 ) 9 à, Co Ed J 5 = J ps 
Po dz Ro dz dt 
(15) 


Admitindo uma solução do tipo: 


O mei ( de) (16) 
: E 
tem-se: 
2% 
e, (17) 
dz t 
99* z 
= — — É! 18 
dt f? um 


Substituindo estas expressões em (15), dividindo 
a expressão obtida por = f” (a solução f'=0 


não interessa, por corresponder a espessura 
constante, que é a solução equivalente a t +» 9), 
virá : 


(19) 
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cujas soluções são: 


T 
* o | 
(é) * — ee nt 


— 2ogcos: 


TT Vi Mo a 
«Vaio rt em 
| Z2pgcosÊ, ogcosb t 
das quais só a solução positiva interessa. 

Finalmente, a espessura à segue a seguinte lei; 


(aa 
29 g cos 
= ; . E = 
+ V (sm) depara ge na A) 
Z2Pgecosb /  B cos É t 
expressão que satisfaz às condições limites : 
À — z=0 
N (a 
q = — — — (22) 
cg cos) 
BD t » 09 
p=. *o (23) 
pg cos |) 
C— Caso particular do fluido newtoniano : 
[2,=0 
(24) 
aum 
Neste caso: 
ou 
E (25) 


nie a 

V pgcosó + 

Coincide com expressão apresentada em [1). 
3. MODELO DE OSTWALD-DE WAELE 


A lei é: 


- (26) 


Como o gradiente é positivo, a equação (26) 
pode exprimir-se : 


Cuz =-—m (E ) 
dx 
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(26") 


3.1. Balanço material e expressão da velo- 
cidade média 


A realização de um balanço material sobre 
a porção tracejada (fig. 2) fornece, depois de 
divisão por AZ At e passagem ao limite [1): 


| R dj 
L(cv>)tÃ=o (27) 
dz ot 
x 
NIVEL INICIAL DE LIQUIDO / 
z ç + 
Z+AZ 


(ja 


NIVEL AO FIM DO 
TEMPO t 


Fig. 2 — Perfil do escoamento, numa secção 
normal à placa, 


De (3) e (26') resulta, explicitando em ordem 


dv 


a cl , integrando e obrigando à condição 


fronteira de ser v, = O para x =0: 


1 na+1 
v, = (*EcosÉ) a 2; [== 
: : cd n+1 l 
Bra. 
à qual corresponde um valor médio : 
; 1 n+1 
“pg cos L É 
Ear ss PRESS ) Do (29) 
2n+1 m 
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3.2. Determinação do valor de à 


Substituindo (29) em (27), admitindo uma 
solução do tipo: 


o (zt) = W (z) o (t) (30) 


e separando variáveis, obtém-se: 


1 1 
eg cosfin yjndv = 1 df c 
mm d> 2n+1 dt 
n 


(31) 


equações que integradas e sujeitas às condições 
aos limites (espessura nula no ponto z=-0, para 
todos os instantes além do inicial e espessura 
infinita no instante inicial) dão origem a: 


n 


1 
W = Kg n+1 Egas Jr | do (32) 
n Pg cos, ” 


e 
EE. 
5 =| c, "Tê. | n+1 (33) 
À n 
e portanto: 
q n 
og cos É vê 


que satisfaz ao caso particular de fluido newto- 
niano (m=v,n=1): 


o E E 
: = t (25) 
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NOMENCLATURA 


g — aceleração da gravidade 
m,n — parâmetros da lei de (Ostwald-de 
Waele 
v, — velocidade do fluido segundo z 
“v, > — valor médio da velocidade do fluido 
segundo z, nas zonas de gradiente 
variável] 
v — valor da velocidade v, nos pontos 
z G [9*,9] 
é — inclinação da placa em relação à ver- 
tical 
9 — espessura da camada de fluido 
* — espessura da camada de fluido de 
Bingham em que o gradiente de velo- 
cidade é variável 
: — densidade do fluido 
u — viscosidade do fluido 
“., T, — parâmetros da lei de Bingham 
Z., — força de corte que se exerce num plano 
normal a x, segundo z 
x, Y, Z — coordenadas rectangulares 
t — tempo. 
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OPTIMIZAR A UTILIZAÇÃO DUM SISTEMA HIDRÁULICO 


APLICANDO A PROGRAMAÇÃO 


RESUMO 


Neste artigo apresentam-se alguns dos pontos principais 
dum método desenvolvido para melhorar a utilização duma 
bacia hidrográfica cuja finalidade primerdial é tornar 
possivel o abastecimento de água a certa região, Para poder 
definir uma política óptima de regularização dos caudais 
através do controlo dos níveis de água dos reservatórios 
incluídos nessa bacia será necessário atender à possivel 
ocorrência de cheias, por vezes causadoras de grave: pre- 
juízos. 

4 aplicação da Programação Dinâmica a este problema 
torna possível a sua resolução, 

Este estudo foi feito considerando como exemplo a bacia 
do rio Dee, no Pais de Crales, para a qual se apresentam 
alguna dos resultados obtidos, 


1. Os benefícios, os custos e a complexidade 
crescentes dos sistemas hidráulicos têm vindo a 
acentuar a necessidade de optimizar não só o seu 
projecto mas também o modo como são utiliza- 
dos. Em muitos problemas de optimização do 
ambito da Engenharia Civil tem prevalecido a 
tendência de considerar como dados inalteráveis 
os objectivos do empreendimento (isto é, os seus 
benefícios) considerando apenas como variáveis 
os meios utilizados para os atingir e que estão 
associados a certos custos. 

Esta atitude (de sub-optimização) conduz, com 
frequência, a soluções bastante afastadas das 
mais económicas pelo que se tem vindo a acen- 
tuar a necessidade de o objectivo do projecto ser 
também tomado como variável passando a con- 
siderar tal realização como uma parcela dum 
plano mais vasto e global de empreendi- 
mentos. 

Para contabilizar os custos e os benefícios 
resultantes da interacção estabelecida entre o sis- 
tema e o exterior é essencial, antes do mais, 
definir criteriosamente as suas fronteiras, o que 
por vezes é mais nebuloso do que se poderia 
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DINÂMICA 


por: LUÍS VALADARES TAVARES 
Engenheiro Civil (1.S.T.) 


Master in Operational Research 
(M, A., LANCASTER) 


SYNOPSIS 


In this paper, it is presented a method to improve the 
utilization of water resource systems having as main purpose 
to supply water. To determine the optimal policy for reser- 
coir operations, the occurrence of foods has also to be 
considered because their effects are sometimes very important. 

Dynamic Programming is successfully usel to solve 
this problem. 

Using data coming from the river Dee basin, in North 
Wales, some of the derived results are herein presented, 


julgar. É clássico o exemplo de pretender avaliar 
os benefícios provenientes de certa obra de irri- 
gação para o que se consideraram não só os seus 
próprios benefícios mas também os provenientes 
de outros sistemas associados. 

Determinar o valor dos custos (C) e muito 
especialmente o dos benefícios (B) é tarefa bas- 
tante difícil até porque estes assumem geralmente 
formas muito diversas. Só se uns e outros forem 
exprimíveis nas mesmas unidades poder-se-á 
definir uma função utilidade do sistema (U == B — 
— C) pois no caso contrário será apenas possível 
definir a sua eficácia (E = B/C). 

Os estudos relativos à utilização dos sistemas 
hidráulicos embora se apresentem sob formas 
muito diversas têm sempre um mesmo objectivo: 
estabelecer sistemas de regras que permitam ao 
longo do tempo tomar as decisões mais aconse- 
lháveis em função da informação disponível. 
Verifica-se geralmente que: 


a) A previsão dos efeitos das várias decisões 


possíveis deve ser feita atendendo à esto- 
casticidade do sistema; 
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Tal como é habitual nos estudos de pla- 
neamento será necessário definir um limite 
no tempo, horizonte até ao qual se façam 
essas previsões ; 

A informação disponível varia com o ins- 
tante considerado ; 

A utilização do sistema implica tomadas 
de decisão sequenciais ; 

Geralmente apenas é possível definir a efi- 
cácia do sistema e não a sua utilidade. 


b) 


c) 


d) 


e) 


Em muitos dos projectos realizados neste do- 
mínio tem sido essencial a directa colaboração 


-BAÍA DE LIVERPOOL ——— 


a o 


po 

dani nm 

— CELYN 
( = 

BALA 
o e «di nes ud 
paço 
ad: OP 
Fig. 


entre técnicos de Investigação Operacional e de 
Hidráulica (!). 


(1) Como exemplos são de destacar as investiga- 
ções feitas na Califórnia pelo Stanford Research Institute 
e pelo Water Resources Centre (University of Califórnia) 
e em Inglaterra pelo Departament of Operational Research 
(University of Lancaster), Este trabalho procura apresentar 
alguns dos pontos principais dum projecto realizado neste 
último departamento por J. C. Wilkinson (orientador), 
D. J. Ashton, o autor deste artigo, A. R. Waters e R. 
Williams. 
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2. O sistema em estudo é constituído pela 
bacia hidrográfica do rio Dee. Este rio nasce nas 
montanhas de Merionethshire no País de Gales, 
possui uma bacia hidrográfica com cerca de 400 
milhas quadradas de área e vai desaguar na Baia 
de Liverpool (Fig. 1). Três importantes reserva- 
tórios contribuem com as suas águas para os 
caudais escoados no Dee: Bala, Alwen e Celyn. 
Os dois primeiros são naturais enquanto que 
o terceiro é artificial. Só Bala e Celyn (com as 
capacidades de 15 e 81 milhões de m/, respecti- 
vamente) são controláveis. 


— aq 
Ea ba 
- lr 
Ta = 
NS 
CHESTER, 
Aw tomada 
diva 
DE e 
Em a Fal 
DO A 
A PE ”, 


v I MILHA 


l 


Os caudais escoados na secção A deste rio são 
decomponíveis em duas parcelas: 


E=E, + E, 


E — Escoamento total 

E, — Escoamento proveniente das bacias pró- 
prias de Bala e Celyn 

E, — Escoamento proveniente da restante área 
da bacia hidrográfica do Dee, a mon- 
tante de A. 
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Enquanto que o valor de E, é regulável o de 
E, não o é. 

As formas por que este sistema hidráulico 
é útil são de natureza muito diversa, tal como 
é aliás usual. A finalidade principal consiste em 
garantir o abastecimento de água a certa popu- 
lação e os benefícios secundários são numerosos 
(constituir um importante meio de transporte, 
fomentar a indústria do turismo na região, etc.). 
Os principais benefícios negativos são os pre- 
juízos provocados pelas inundações. Estes prejuí- 
zos englobam perdas materiais e efeitos de 
carácter sócio-político., Com efeito, esta região 
do País de Gales é das que apresenta aspirações 
mais vivas de emancipação regional e como as 
autoridades responsáveis pelo controlo dos níveis 
de água nestes reservatórios são inglesas, a ocor- 
rência duma cheia é sempre motivo para a popu- 
lação local apresentar inúmeros protestos, por 
vezes de forma violenta. 

Diminuir a frequência e a dimensão das inun- 
dações exigirá que se mantenham nos reserva- 
tórios os níveis de água tão baixos quanto pos- 
sivel. Em contrapartida, assegurar maiores cau- 
dais para abastecimento implicará manter em 
cotas máximas esses níveis. Pretende-se, portanto, 
obter um compromisso óptimo entre dois objecti- 
vos em conflito. 

A análise hidrológica da bacia mostrou que 
as cheias podem ser provocadas por períodos de 
precipitação intensa de pouco mais de quatro 
horas enquanto que as secas exigem períodos 
de estiagem de dois ou três meses. 

Atendendo à existência desta diferença entre 
escalas no tempo para cheias e secas e ao facto 
de as autoridades locais entenderem que o abas- 
tecimento de água é um objectivo manisfesta- 
mento prioritário, adoptou-se a seguinte atitude; 


1) Maximizar o caudal de abastecimento, for. 
necido com uma certa garantia, estabele- 
cendo quais os níveis mínimos a manter 
nos reservatórios ao longo do tempo, 
Obter-se-á, deste modo, uma política de 
regularização a médio e longo prazo. 

2) Respeitando os níveis limites anteriormente 
determinados, definir uma estratégia de con- 
trolo a curto prazo do sistema que mini- 
mize a ocorrência de cheias. 


Neste artigo apenas se apresenta a resolução 
do primeiro ponto. 
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3. O que há de controlável neste sistema hi- 
dráulico é o modo de regular o nível de água em 
Bala e Celyn. As suas bacias hidrográficas têm 
características bastantes semelhantes no que res- 
peita à pluviosidade e as suas áreas estão na razão 
de 1 para 3. À zona de Bala é (ao contrário da de 
Celyn) uma das áreas em que as inundações têm 
efeitos mais desastrosos pois a cidade de Bala 
situa-se junto a esse reservatório e não dispõe 
de qualquer protecção. A capacidade do reserva- 
tório de Bala representa apenas 15 “jo do total 
armazenável em Bala e Celyn. Bala terá uma 
função de desprezar enquanto reservatório que 
acumule volumes de água significativos para 
o abastecimento de água. 

Tem-se ainda que o problema das secas, como 
já foi dito, se situa numa perspectiva de médio 
e longo prazo ao contrário do das cheias. Por 
todos estes factos, é recomendável que se con- 
sidere Bala como um reservatório destinado ape- 
nas à regularização de ondas de cheia durante 
períodos de tempo (curtos) em que se verifiquem 
escoamentos máximos. No estudo da minimização 
da frequência de ocorrência das secas, ignorar 
o reservatório de Bala traduzir-se-á na introdu- 
ção dum erro desprezável que será sempre 
a favor da segurança. Para este efeito, o sis- 
tema inicial será redutível ao da figura 2 
E, = E cE,= E — Ec). 


Ec 


Fig. 2 


Os dados de que se pode dispor são escassos. 
Conhecem-se os escoamentos naturais mensais 
em Bala e Erbistock (secção A) durante 11 anos 
e em Celyn durante 2 anos compreendidos entre 
os primeiros, 

Os dados referentes à pluviosidade apenas 
permitem concluir que o período dos 11 anos 
referidos tem características médias. 

Verificou-se que, a menos de desvios extrema- 
mente reduzidos, existe a seguinte relação : 


Ec = 0,25 E, 
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Sendo: E, — Escoamento mensal em Bala 
E — Escoamento mensal em Celyn 


Consequentemente, o volume de água acumu- 
lado durante um mês em Celyn será: 


V=0,25E,—D, 


Sendo: D, — Pedido mensal de água armaze- 
nada, para abastecimento. 
D — Pedido mensal de água, para abas- 
tecimento. 


Análise estatística dos dados permitiu concluir 
que In E, segue uma distribuição normal com 
média igual a In (k E, + B) sendo k, Be a variân- 
cia constantes. 

Como era previsível, o valor do pedido de 
água armazenada depende, portanto, do escoa- 
mento proveniente da zona da bacia hidrográfica 
de águas não regularizáveis. 

Todavia, os sucessivos valores mensais de 
E, não são independentes pelo que não será lícito 
considerá-lo caracterizado apenas por uma lei de 
distribuição estatística. O modelo apropiado é a 
consideração de E, como um processo estocás- 
tico realizado num espaço de dimensão igual a 12, 
em que os valores de E, ao longo dos 12 meses 
constituem uma família de 12 variáveis aleatórias 
definidoras do processo. Foram assim estimadas 
12 distribuições para os escoamentos verificados 
nos 12 conjuntos de 11 meses. O tipo de distri- 
buição que se mostrou mais adequado foi o da 
distribuição Gama com parâmetros variáveis de 
mês para mês, 

Utilizando os escoamentos não regularizáveis 
dessazonalizados obtiveram-se os seguintes coefi- 
cientes de correlação (Tabela 1). 

Usando estes novos valores é possível afirmar 
que E, estará distribuído segundo uma certa lei 
com média (u ) e variância (02) constantes e que 
a distribuição de V será condicionada pelos valo- 
res de E,. Aplicando a fórmula de Bayes foi 
possível estimar a distribuição de V com pará- 
metros uy e oy” constantes. Por definição, a seca 
resulta da incapacidade de satisfazer certo pedido, 
isto é, de o pedido de água armazenada não 
poder ser satisfeito o que se verifica quando o 
nível de água no reservatório se encontra abaixo 
de certo valor mínimo. Pode-se portanto definir 
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TABELA 1 


Função de autocorrelação de escoamentos mensais 
em Bala (1957-1967) 
(Meses) Usando valores Usando valores 
reais dessazonalizados 
() 1.000 Odo 
| 0.273 O.010 
2 (0.003 — 0,007 
3 RNUE 0,004 
| — (0,107 — (0,003 
É — (1,314 — (1018 
() - 0,201 U4S 
7 — (1,210 005] 
3 - 0,009 QUO 
Uy — UOUS — D.105 
ly — OMNI — (100 
1 (1,385 UIT 
12 (1,373 — 0015 
13 (1.240 — (0,039 
14 (1,250) (1,248 
15 — 043 -— QUSU 
Lt — U.I54 — (1,113 
17 | — 0,202 0.089 
E | — (0.336 — 0.013 


seca como a descida do nível de água existente 
no reservatório abaixo duma certa cota limite. 

4. O objectivo deste trabalho é, pois, o 
seguinte : 

Para um certo nível de água B, valor de pro- 
babilidade p, e prazo h, determinar uma série 
de níveis de retenção ao longo do tempo, L (t), 
tais que, para qualquer instante t, em que o pre- 
sente nível de água x (t) seja L (t), a proba- 
bilidade de que o nível de água permaneça supe- 
rior a B, pelo menos até t + h, seja igual a p. 

A resolução deste problema será conseguida 
atravez da Programação Dinâmica que, como se 
sabe, é baseada no Princípio de Bellman. 

Discretizando, consideram-se 12 intervalos de 
tempo coincidentes com os 12 meses do ano. 

Pretende-se L, (i=1,....,12) tais que, 
se x, = L, no mês j, então a probabilidade 
de x, > B, para todo o K inferior ou igual a 
| + h superior a j, será igual a p. 

Seja f(x,i,n) a probabilidade de que não 
ocorra seca nos próximos n meses, sendo x 
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o nível de água no mês i (1,....,12). Sendo y, 
o volume de água útil que entrou no sistema no 
mês i, 6, (y) a sua função densidade de 
probabilidade e P, (y,) a sua função de distribui- 
ção, existirá a seguinte relação de recorrência: 


Ljry—X 
f(x,i,n)= / f(x+ypi+In—1)6,(y;) dy + 


Pr, — (O) 
+f(L, yo n+7,n—1) 1-9, (L,,4—>)] 
(sei=12,i4+1=1) 


Sendo Z,=y,—S, em que S, representa o valor 
médio de y, 


L, 4 |—x—S, 
Eb tempo | fix 5, + Zul + 


JU — CO 
+ 1,n—l1) d(Z,) dz, + f(L, ppi + Tni=—I)>< 
se [1 + dp (L,,,/—x—59] 


Começando com f(x,i,0) = 1 podem-se calcular 
os valores até f (x,ij,n) para cada x (B< Xx< Lj), 
para cada i (i=1,....,12) e para valores 
sucessivos de n (n==1,...., N). 

Calculam-se assim as probabilidades de não 
haver seca em função dum conjunto de níveis 
de água no reservatório. O que se pretende é 
obter os níveis correspondentes a uma certa 
probabilidade de seca. A relação p, = f (L, i, N) 
não permite ser analiticamente invertida para 
obter L,=g (p,i,N) 


[L, | [Pi | 
Sendo r — É ep = p, 
(Lo) [Pio] 


procedeu-se do seguinte modo : . 
Determinou-se r “ sendo conhecido p e T 
tal que 


pf=T(r*) 


Ts é a solução inicial: p, — T(r,) 


:, será o vector variação de r do 
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Variando i de 1 a 12 pode-se construir a 
matriz D em que 


dp 
o(i) 
A 
Fo(ii 
Então 
12 dp 
k cias à o(i) 
É Pop — > à Tor) Edi 
j=1 0 Toti) 


Portanto p, = D3r, ou r=D"'ir,. 

Substitui-se a solução inicial r, por 1; =, + 
+ dr, 

O processo será dado por terminado quando 
se considerar ser Pk suficientemente aproximado 
de p. 

4. Fazendo B = 897 pés (corresponde a um 
volume de água armazenada de 18 milhões 
de mº), h=-5 anos e p==0.825 obtiveram-se 
os seguintes níveis que contrastam com os que 
tradicionalmente se usaram (fig. 3). 

Usando uma certa série de níveis foram 
calculadas as probabilidades de não haver seca 
em 1, 2, 3 e 4 anos e constatou-se que a curva 
obtida era a que resultava de se admitir que o 
intervalo de tempo entre secas seguia uma lei 
exponencial negativa. Impondo que essa lei tome, 
exactamente, o valor verificado para o período 
de 1,2, 3 ou 4 anos verificou-se que a curva se 
adapta tanto melhor quanto maior é esse período 
(tabela 2). 

Com os níveis óptimos obtidos para p == 0.825 
em 5 anos, foram calculados os valores de p para 
períodos de duração até aos 25 anos. Uma vez 
mais o intervalo de tempo entre secas (AT) 
seguiu uma lei aproximadamente igual à da 
distribuição exponencial negativa (fig. 4). 
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NÍVEIS LIMITES EM CELYN 
(expressos em pés) 


— — — — —— — 


970 


960 


950: 


| a mia Gê —p Valores usados 
| tradicionalmente 


—— so Valores propostos 
940 ' | ! ! | 


Vig. 3 
TABELA 2 


Probabilidade de não ocorrência de seca utilizando o método apresentado e a lei de distribuição exponencial negativa, 
Apresentam-se, como exemplos, os valores obtidos para o mês de Janeiro e Junho. 


JANEIRO 
| | | 
Tempo (Anos) | | | 2 3 4 | 
= | | 

Método apredentado cce mngre ses em 0,981] | 0.937 (1,802 (1,848 

| 
Exponencial negativa (t= 1) . cs cuco + +| (0081) 0,062 (1.944 0,926 
Exponencial negativa (t= 2: . ES | 0,008 (0,037) O.9U8 1,878 
Exponencial negativa (t= 3... Es a a | 1,003 (1,027 (11,892) (1,850 


Exponencial negativa (t = 4) ER E AS A | 0.960 (1.921 (1.884 (0,848) 


— —— ————— o o e a —— a SS e e tt TT 


JULHO 
| 
Tempo (Anos) | | 3 

| 
Método apresentado . . .. SERES ÉS (OU | NOTOU (LD (SSO 
Exponencial negativa (1 = 1) cce cll cc. (1,000) (U.908 (1.997 RT 
Exponencial negmtiva (=2) scenes ces (085 (0.970) (1955 (041 
Exponencial negativa (t=3) +. ccccscsc..s, “4 1.074 0040 (0,025) 0.901 


Exponencial negativa (( = 4) .. cc cc cc cc. U9QU | 1,938 | 0.000 (0,880) 
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Fig. 4 
Pode-se então concluir que o número de secas 2 


verificadas num pesíodo de tempo T seguirá a 
lei de Poisson com média u ==)T. Sendo T 


igual ao valor médio do intervalo de tempo entre ” 
secas (“AT), a probabilidade de não ocorrer 
nenhuma seca em t anos será 

4. 


st ds 


Só devido a esta conclusão existirá uma corres- 
pondência biunivoca entre p (t) e T. 

Como neste caso t=-5 anos e p= 0,825, 9 
os níveis obtidos correspondem a ocorrência de 
secas com um intervalo de tempo médio de 25 
anos. Consequentemente, a probabilidade de se 6. 
verificar pelo menos uma seca em 25 anos será 
de 0.63. 
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Conception, Études et Construction de Parcs Souterrains 
pour voitures à Paris 
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ment et des Travaux Publics) 


Alternative Routes to Steel 


The Influence of Microstructure on Low-Temperature 
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Meeting of the ASTM Committee B-8 on Electro-Depo- 
sited Metal Coatings and Related Finishes 
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garo com traduções simultâneas 
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ESTUDO DO FLUXO EM MEIOS ANISÓTROPOS PELO 
METODO DOS ELEMENTOS FINITOS 


RESUMO 


Aplica-se o Método dos Elementos Finitos do cálculo das 
forças de percolação em meios constituídos por camadas 
estratificadas, com fronteiras seguindo uma geometria 
qualquer e direcções de permeabilidade variáveis, Apresen- 
tam-se alguns exemplos ilustrativos e o cálculo da linha de 
saturação numa barragem de terra, 


1— INTRODUÇÃO 


Em meios isótropos, homogéneos e de frontei- 
ras particularmente definidas, o cálculo da per- 
colação, obtido pela solução analítica da equação 
diferencial de Laplace ou pelo método gráfico 
baseado nas propriedades desta equação, é per- 
feitamente eficaz e correcto. No entanto, a não 
ser que se admita uma idealização bastante gros- 
seira, o solo, na maioria dos problemas da prática, 
jamais assim se apresenta, nem mesmo nos locais 
onde certas zonas são elaboradas pela mão do 
homem, como nas barragens de terra. 

Houve, portanto, necessidade de procurar méto- 
dos de cálculo mais gerais. 

O desenvolvimento e a crescente utilização dos 
computadores tornaram possível a formulação de 
novas técnicas para solução rápida e fácil de 
problemas de difícil e mesmo impossivel solução 
analítica. 

Um passo à frente foi dado com a utilização 
de métodos numéricos. Inicialmente, aplicou-se o 
método das diferenças finitas. Mais recentemente, 
com o desenvolvimento do método dos elemen- 
tos finitos, verificou-se sua perfeita adequabili- 
dade à solução de problemas deste tipo, possibi- 
litando o cálculo da percolação, dos casos mais 
simples aos mais complexos, sem praticamente 
nenhuma diferença no tratamento teórico ou no 
tempo computacional. Assim, várias zonas de 
diferentes solos, heterogéneas, anisótropas e com 
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SYNOPSIS 


The Finite Elements Method is applied to the compu- 
tation of seepage forces in atratified media with any 
boundaries and directions of permeabiliiy. Some cases are 
presented such as, for instance, the calculation of the 
phreatic line of an earth dam. 


total arbitrariedade na geometria de suas fron- 
teiras, são tratadas segundo os mesmos princi- 
pios, e, conforme adiante se verá, os resultados 
obtidos são tão aproximados quanto se julgar 
necessário. 

Pimeiramente, considere-se o conceito, introdu- 
zido por Ferrandon, de matriz de permeabilidade. 
Estabeleceu ele uma matriz de coeficientes de 
3x 3, simétrica, e admitiu que em qualquer mate- 
rial existem 3 eixos principais na direcção dos 
quais há uma colinearidade entre a percolação e 
o gradiente de pressão. Os coeficientes de per- 
meabilidade, obtidos para um volume de solo 
particular segundo a direcção destes eixos, des- 
crevem a permeabilidade neste volume. 

Estas direcções, por serem determináveis em 
laboratório, são consideradas conhecidas. 

A partir daí, a lei de Darcy, bem conhecida e 
estudada, pode ser generalizada para a percola- 
ção em meios anisótropos e heterogéneos. 

Os métodos gráficos e analíticos são clara- 
mente impraticáveis quando na zona percolada 
existem várias direcções principais de permea- 
bilidade. 

Todo o estudo aqui desenvolvido visa os casos 
da prática em que é procurada a solução de pro- 
blemas a duas dimensões, sendo simples e ime- 
diata a generalização para as três dimensões, se 
se desejar. 

Além disto, materiais como rochas e solos são 
geralmente estratificados, existindo igualdade das 
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características fisicas segundo o plano de estra- 
tificação. Na fig. 1, os eixos x e z são paralelos 
ao plano de estratificação. 


Y 


2 — FLUXO EM MEIO ANISÓTROPO 


Designa-se por v o vector da velocidade de 
percolação, com componentes segundo os eixos 


x e y, respectivamente v, e v,, 


e por grad H o gradiente hidráulico da carga 
total H, 


dH 
Ox 
dH 

0y 


grad H = 


A relação linear mais geral entre estes dois 
vectores é estabelecida pela lei de Darcy: 


(v) = — |K| (grad H) (1) 
onde K é uma matriz da forma 


K | as Ra Ko 
K K 


yvx vy 


obtendo-se, pelo príncípio de conservação da 
energia, 


Considerando novos eixos x”, y', fazendo um 
ângulo &) com x, y, o vector velocidade segundo 
estes novos eixos é definido como 


238 


(v] = |R| (vi (2) 


onde R é a matriz de rotação, cujos termos são 
os co-senos directores dos ângulos que entre si 
formam x, y e x”, y”, ou seja: 


pa (3 sen O 
— sen O cos €) 
Da mesma forma 
fgrad H') = |R| [grad H) (3) 


Combinando (1), (2) e (3) entre si, obtém-se 
(v)=—I|R| |K| |R7! fgrad H) (4) 
E daí, como 
IK'| fgrad H) = — |R| |K| RT! (grad H) 
virá 
K'|=— RI |K| RT 


E é sempre possível determinar um ângulo O 
para o qual K' se reduz a uma matriz de diagonal, 
onde : 


Í ! 
Rai 


0 Ke 


a 


(grad H') 


As direcções dos eixos segundo os quais a 
matriz de permeabilidade se resume a uma matriz 
de diagonal definem as direcções dos eixos 
principais de permeabilidade. 

Num material estratificado, tal como rocha ou 
solo, os eixos principais estão localizados segundo 
o plano de estratificação e normalmente a este. 

A «equação de continuidade» do escoamento 
de um líquido incompressível através de um 
material poroso, que não sofre variação de volume, 
é bem conhecida: 


Desde que as direcções x e y coincidam com 
as direcções principais de permeabilidade, esta 
equação pode escrever-se sob a forma 
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denominada «equação de fluxo». 
3 — SOLUÇÃO DA EQUAÇÃO DE FLUXO 


Conhecida a carga H actuando em determinada 
região, estabelecidas as fronteiras de escoamento 
e características de permeabilidade das diversas 
zonas, há apenas uma forma de escoamento, e as 
pressões são bem definidas em cada ponto da 
zona percolada. Existe portanto apenas uma 
solução para a equação do fluxo. 

De acordo com os princípios do cálculo varia- 
cional, a solução desta equação é equivalente 
à determinação da função H que minimiza o 
seguinte funcional, tomado sobre toda a região 
do fluxo: 


, + 9 
2 7 dx) 
= 
Te Rss = | axdy (5) 


O método dos elementos finitos fornece uma 
técnica para determinar os valores de H em 
alguns pontos do domínio que tornam mínimo 
o referido funcional. 

Não é objectivo deste trabalho mostrar deta- 
lhadamente o que é o método dos elementos 
finitos, pois que sobre isto já existe extensa 
bibliografia. 

Faz-se no entanto um breve relato da aplicação 
do método a este problema. 

A região do fluxo é dividida num conjunto de 
elementos finitos triangulares (Fig. 2). 


cade 


TECNICA Nº 405 


Considera-se um elemento típico de nós i,j,m, 
e, em toda a zona limitada pelo elemento, 
supõe-se ser H contínuo e função apenas das 
coordenadas (x,y) do ponto. 

A função mais simples é definida pela relação 
linear : 


H=A+Bx+Cy (6) 


onde A, Be C são constantes. 

Pode-se agora, em função das coordenadas dos 
nós i, jm, (X, Y), (xy Yj), (Xmr Ym) Tespectiva- 
mente, estabelecer um sistema de 3 equações 
a 3 incógnitas e eliminar A, B, e €, obtendo H 
como função das coordenadas dos nós e dos 
HH nestes nós. 


H=A+Bx, + Cy, 
H, = A + Bx, + Cy, 
H, = A+ Bx, + Cym 


em notação matricial, 


H, 1x, Y; H 
H [=|1x, y |jB 
Ela 1Xm Ym Lc] 


Por intermédio destas equações calculam-se 
os valores A, B, e C, os quais, substituídos 
em (6), permitem obter H: 


H=|(a+bx+cy) (a +b;x+cy), 
HI 
H 


(a, +bax + cy) 1H | (7) 
E 


m 


sendo a,, b, e c; definidos por: 


Ee 
| 
2 a 
ea 
| 
| 
as 
E cd 


Nestas expressões À é a área do triângulo, e 
os coeficientes, a, b,, c e am Das Cu SÃO 
obtidos variando i, j, m, ciclicamente no sentido 
anti-horário. 

Os valores (ainda desconhecidos) de H nos 
nós do elemento definem a função completa e 


unicamente através da região percolada. 
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Derivando o funcional F em relação aos HH 
nodais e igualando a zero, resulta um sistema 
de equações lineares. 

Portanto, derivando F no elemento (F") em 
relação a H,, H e My vem: 


o ES ue 9H J “ES 
me K Eni 
9H, ff = dx dH VET ii 


d H E) dH 
Um Ksy dy d H, d y d, d, 


Substituindo H pela sua expressão dada em (7), 


H, 
E ff x Ibi by b Im IN + 


DH, do | KN m | | 
Em 
H, 
+ Ko] Er 0746] 1H || |d,d, 
H,, 
Como 
aci Rand b,, b,, bme Cs Co CH HH, 


são independentes de x e y, tem-se 


9 ES 
5H. as Rs b;b, + K,y€; C,, K ex D; b, + 
i 
cmi 
te Roo 8 Cio Sax O D+ Km G4|><4H,+ À 
EN 


Da mesma forma, obtém-se a derivada de EF 
em relação a H e H,,o que dá: 


- simétricos 
| d Hº | d H, 
9 Fº 
d m 
| 
Xi H |=8X5S,,><H 
EM 
z40 


A soma das contribuições obtidas pela minimi- 
zação de cada nó, em função dos HH que definem 
a continuidade em cada um dos elementos que 
têm em comum este nó, resulta uma equação 
linear. 

Cada equação é da forma 


2H 2/4 in se 5], =0 


H, representa o H no nó de um dos | elementos 
os que contêm o nó 1. 


e . . 
S;n representa os termos da linha i de uma 


onde 


das | matrizes semelhantes a S,,, formada pelos 
elementos para os quais o nó | é comum. 

À in representa a área de um destes | elementos. 

Forma-se assim um sistema de equações em 
que as incógnitas HH são tantas quanto são 
os nós. A simples observação mostra-o como 
sendo simétrico, em banda, e de segundo mem- 
bro nulo. Introduzindo se neste sistema os HH 
conhecidos nos nós de fronteira, e resolvendo-o, 
obtêm-se os HH em todos os nós dos elementos 
em que foi dividido o meio percolado. 

Observa-se ainda que quanto mais estreita for 
a malha de elementos finitos mais precisa será 
a solução. 

O processo descrito é semelhante ao do cál- 
culo dos esforços em estruturas pelo mesmo 
método, e o computador é utilizado segundo 
o mesmo roteiro. 

Na fig. 3, mostra-se uma aplicação do método, 
para um caso bastante particular e simples. 
A fig. 4 apresenta a solução para um caso bas- 
tante geral, notando-se claramente nas frontei- 


-——— a mt, a mm am eme me mma Re 


Ka Did; + Ko CG Go Kas Dj dy + Kss & Go Ku ba d+ K,50€ 


K x D,D; + Kg € Cir Kyx Dm D; + Ky, Cm 6; | >< 


Ka b,, Db, E K sy € | 


EH Cm | 
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DE cd 


NOTA -AS LINHAS TONTINHAS MEPRESENTAM 
EQUIPOTENCIAIS CUIA CARGA ÉSIA 
INDICADA 


as ão 


a 


Fig 3 - Equipotenciais 


ras das zonas a «refracção» da equipotencial 
devido à mudança de permeabilidade. Nota-se 
também que apenas a existência de 2 zonas dis- 
tintas não justificaria a perturbação na forma 
das equipotenciais. É, sem dúvida, nas de carga 
4 e 5 que mais fortemente se evidencia a influên- 
cia conjunta dos vários factores intervenientes: 
duas camadas distintas, não isótropas, frontei- 
ras não lineares e direcções principais de per- 
meabilidade variáveis. 

Quando um dos limites de percolação não 
é conhecido à priori, embora fisicamente bem 
definido, como é o caso da linha de saturação 
nas barragens de terra, torna-se necessário veri- 
ficar se o campo estabelecido para o fluxo pelo 
lançamento arbitrário da malha de elementos 
finitos triangulares está compreendido e limitado 
pelos limites da região de percolação. 

O procedimento adequado é o seguinte: 


a) Lançada uma rede de elementos finitos, cal- 
culam-se os HH nos pontos nodais, e determi- 
na-se uma linha de saturação. 
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b) Verifica-se se os elementos finitos adopta- 
dos estão geomêtricamente coincidindo com 
os limites impostos pela linha de saturação cal- 
culada. 

c) Caso não estejam é necessário refazê-los 
na zona próxima à linha de saturação, aproxi- 
mando os limites da malha aos limites da linha 
de saturação anteriormente calculada. 

d) Isto deve ser repetido até que haja uma 
coincidência adequada entre os limites da malha 
e a região de fluxo. 


Normalmente, 3 ou 4 iterações são suficientes 
para se obter a necessária precisão, conforme 
se vê na fig. 5. 


4 — CONCLUSÕES 


O método descrito permite uma aplicação prá- 
tica a grande variedade de problemas. Sem neces- 
sitar de programas especiais, além daqueles utili- 
zados no cálculo de estruturas, mostra grande 
eficiência no que diz respeito a tempo computa- 
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cional. Além disto, permite ainda uma extensão 
a problemas tridimensionais utilizando-se tetrae- 
dros em vez de triângulos. 
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4 HERÓIS DA LUBRIFICAÇÃO > 


Punctilius 


No ano 56 a. €., o imperador Júlio César, que há muito 
pensava escrever um livro sobre a conquista das Gálias, 
conquistou-as mesmo e escreveu o livro. Se Você estudou 
Latim, bem sabe (olá se sabe!) a que livro nos referimos 

Mas, como muitos dos grandes desse mundo de antes 
de Cristo, César esqueceu-se de mencionar na obra o papel 
desempenhado por Punctilius — chefe de manutenção dos 
exércitos romanos que, afinal, tornou possivel tal vitória e, 
por conseguinte, o célebre livro. 

No fim da semana que precedeu a invasão, a incle- 
mência dos elementos (sempre mau tempo nos fins de se- 
mana!) deitara por água abaixo a eficácia das legiões 


Punctilius 


romanas. A humidade enferrujara os metais. O sol (deve 
ter sido no sábado, porque não há sábado sem sol) viera 
depois secar os apetrechos de couro. As lanças, as espadas 
e os arietes tinham os gumes rombos. As rodas dos carros 
de guerra chiavam como carros de bois minhotos. 

Os problemas de lubrificação que afligiam as legiões 
romanas eram enormes. César viu-se, bem depressa, enter- 
rado em papiros que lhe exigiam soluções imediatas. 

De cabeça perdida, como qualquer mortal, César 
mandou chamar o seu chefe de manutenção — «Punctilius», 
disse, «vai comprar lubrificantes melhores, senão a gente 
nunca mais passa o Rubicão». 

Punctilius, contudo, chamou, com a devida vénia, a 
imperial atenção para o facto do problema não se resol- 
ver só com lubrificantes. Seria necessário, também, organi- 
zar e dirigir a sua utilização por todo o imenso e disperso 
exército de invasão. 

«Só a MOBIL nos pode ajudar», disse Punctilius. 

«Eles possuem um sistema de processamento de dados 
designado por MI/DAC — que é ideal para a lubrificação 
e manutenção em operações complexas». 

Com a colaboração de Punctilius, os técnicos da MOBIL 
criaram um sistema MI/DAC para o exército romano. 

Os dados relativos às necessidades de lubrificação 
eram imediatamente transmitidos à Secção de Ábacos de 
César, que logo imprimia mapas com recomendações e ren- 
dimentos de lubrificantes. E, ainda por cima, fornecia um 
pergaminho de controle sobre o custo e a eficiência do 
programa inteiro. 

Com métodos de processamento automático, eficiente- 
mente executados, o exército romano, bem lubrificado, 
derrubou a resistência da Gália num ápice. E César escre- 
veu o livro. 

Os chefes de manutenção dos nossos dias podem tam- 
bém, amiúdo, defrontar problemas semelhantes. 

Não basta ter lubrificantes bons. 

Pode ser necessário, igualmente um programa de lubri- 
ficação que abranja todos os casos, um programa capaz 
de resolver todas as exigências de uma produção fabril em 
larga escala. 

Os chefes de manutenção acabarão por verificar como 
um sistema de processamento automático de dados poderá 
resolver eficientemente os seus problemas. 

Não desconhecem, também, que o sistema MI/DAC 
é um dos muitos que a MOBIL criou e pôs à disposição 
dos seus clientes para que eles possam obter o máximo 
rendimento dos produtos que ela lhes vende. 

Punctilius tinha razão ao dizer: 

«A César o que é de César, mas exijam MOBIL!». 


Mobil 
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PAIVA LVA LARA 
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AÇOS 


C. D.T. 669,14.018.2 


INOXIDÁVEIS MARTENSÍTICOS E FERRÍTICOS 


COMPORTAMENTO E ADIÇÕES ESPECIAIS (*) 


RESUMO 


Veste artigo fas-se uma revisão das propriedades 
mecânicas, reacção aos tratamentos térmicos e precipita- 
ção de carbonetos nos aços com 12 a 15 º/y Cr. Dos aços 
com teor em crômio superior a 15 "/o focam-se os proble- 
mas relacionados com os processos de fragilização a que 
estão sujeitos, e ainda as suas propriedades mecáânicas, 
No que respeita a elementos de liga, analisa-se o efeito 
dos elementos mais efectivos, nomeadamente no aspecto 
de propriedades mecânicas a altas e baixas tempe- 
raturas, 


1. GENERALIDADES 


Os aços inoxidáveis, para além da sua boa 
resistência à corrosão, comportam-se, sob vários 
aspectos, de modo peculiar. Limites relativamente 
bem definidos de composição nos diversos ele- 
mentos determinam comportamentos distintos. 
A reacção aos tratamentos térmicos, a estrutura 
e as propriedades mecânicas são alteradas, com 
incidência, por vezes, na própria resistência 
à corrossão. 

Nesta revisão tenta-se uma síntese dos aspectos 
fundamentais do comportamento destes aços e das 
possibilidades de melhoria de características, com 
base em trabalhos recentes sobre o assunto. Não 
se consideram os aços inoxidáveis austeníticos. 
Esta revisão inclui apenas os aços inoxidáveis 
essencialmente constituídos por ferro, crómio 
e carbono, e estuda-se o efeito dos diversos 
elementos de liga adicionais. 

O elemento principal responsável pela boa 
resistência à oxidação é, como se sabe, o crómio. 
A sua acção é devida à formação de uma fina 
película de óxido de crómio, à sua superfície do 
material. O óxido de ferro só deixa de se formar 


por — M. SOUSA CASTRO 
Eng.º Mec. 1. S. T. 


J. M. ALMEIDA SEQUEIRA 
Eng.º Quím,-Ind. 1,S.T, 


SYNOPSIS 


The mechanical properties, heat treatments and car- 
bide precipitation in 12-15 º/y chromium steels are revie- 
wed. Problems related to the brittleness in steels with 
a chromium content higher than 15 º'y are discussed, 

The effect of the most important alloying elements 
in these steels, with special reference to the mechanical 
properties at high and low temperatures, is also analysed. 


a partir de um determinado conteúdo mínimo 
em crómio. Esse é o mínimo teor que garante 
a inoxidabilidade, e na prática toma-se o valor 
12º/, Cr como limite inferior [1]. 

A fig. 1 representa cortes do diagrama Fe-Cr-C 
e irá servir de base à exposição que se segue. 
A posição das linhas «solvus» e «solidus» não 
está perfeitamente determinada em muitas zonas. 
O crómio retarda considerávelmente todas as 
transformações, o que dificulta a obtenção experi- 
mental das condições de equilíbrio. 

Em função da estrutura, os aços inoxidáveis 
constituídos por ferro, crómio e carbono podem 
ser classificados em três grupos: 


Martensíticos — A alta temperatura são integral- 
mente austeníticos e podem, por 
isso, ser endurecidos por têmpera. 
O teor em crómio varia de 12 
a 18%, e o teor em carbono pode 
atingir 1º. É possível a retenção 
da austenite à temperatura am- 
biente, se o teor em crómio for baixo 


e o teor em carbono alto (ver 2.3). 


(*) Esta revisão faz parte dos trabalhos preparatórios do projecto de investigação iniciado pelo Grupo de Meta- 
lurgia do Laboratório de Física e Engenharia Nucleares sob a designação geral: «Estudo da aplicabilidade de um aço 
inoxidável de origem nacional como material componente de reactores nucleares de potência». 
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e) 


17%Cr 


13 %Cr 


MEERATURA, 


% CARBONO 
Ky las €6 +» Kg 7€s 


Fig. 1 — Cortes do diagrama Fe-Cr-C (Bungardt, 
Kunze e Horn) [6]. 


Ferríticos — A sua estrutura é ferrítica a todas 
as temperaturas. O teor em crómio 
é sempre superior a 16%/. O con- 
teúdo em carbono é de cerca 
de 0,019 para os mais baixos 
teores em crómio e pode atingir 
0,359/ para 30"/ Cr. 

Martensítico — 

— Ferríticos — A sua estrutura a alta temperatura 
é constituída por ferrite e auste- 
nite, em proporções variáveis com 
a composição. 


Quando o conteúdo em crómio ultrapassa 15º/ 
levantam-se problemas específicos, que têm sido, 
aliás, o maior obstáculo ao uso destes aços. 
Far-se-á, por isso, a análise separada dos gru- 
pos 12-15“ Cr e 0 Cr> 15. 


2. AÇOS COM 12-15º/ CRÓMIO 
2.1 — Aços com 0,35“ C 


Com a adição de crómio, o ponto eutectóide 
do sistema Fe-C degenera. Para 12º/o Cr o ponto 
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25 %lr 


que lhe corresponde encontra-se a cerca de 
0,38"/) C e 830 ºC [2]. Segundo outros autores, 
o teor em carbono daquele ponto é de 0,3"/o [3] 
e de 0,35" [4]. 

No aspecto estrutural, um aço com 12º/% Cr 
e cerca de 0,35º/) C é equivalente a um aço de 
carbono de composição eutectóide (0,8"/0 C): 
quando arrefecido lentamente, desde a zona 
austenítica, apresenta uma estrutura integralmente 
perlítica. 

Quando arrefecidos rapidamente desde a zona 
austenítica, estes aços apresentam, à temperatura 
ambiente, uma estrutura martensítica. O carbono 
encontrar-se-á dissolvido na martensite ou não, 
conforme a estrutura existente antes do aqueci- 
mento, e em função do tempo de permanência 
na zona austenítica. Como se disse, o crómio 
retarda todas as transformações. A difusão do 
carbono na austenite é mais lenta nos aços com 
crômio do que nos aços de carbono, e por isso 
é necessário mais tempo para a dissolução dos 
carbonetos. No arrefecimento, esta acção do 
crómio tem como resultado um aumento de 
temperabilidade. Os aços com 12-15"/y Cr são já 
praticamente temperáveis por arrefecimento ao ar. 

Após a têmpera, os revenidos até 500 “C não 
produzem amaciamento sensível (fig. 2). A ligeira 
redução de dureza que se pode verificar até cerca 
de 350 “C segue-se um pequeno endurecimento, 
entre 450 e 500 ºC. Segundo MONYPENNY [3] este 
endurecimento é devido à sobreposição de dois 
fenómenos: transformação, à temperatura de 
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TEMPERATURA DE REVENIDO, C 
Vig. 2 — Eleito da temperatura de têmpera no revenido 


posterior de um aço inoxidável (0,32 /, C, 12,2 */, Cr). 
à - Temperado de goo *C; B -Temperado de 1050ºC [3]. 
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revenido, de austenite intergranular retida em 
martensite, e precipitação fina de carbonetos. 

RONALD e BODSWORTH [5] efectuaram ensaios 
sobre um aço com 11,5º/ Cr e 0,2º/ C e con- 
firmam a precipitação de carbonetos no revenido 
a 500 ºC. A fig. 3 indica a variação das carac- 
terísticas mecânicas com o tempo de revenido, 
segundo os resultados daqueles autores. Encon- 
tram-se igualmente representados os tempos de 
precipitação de dois dos carbonetos de crómio. 
A fase carbonetos será tratada com mais por- 
menor em 4, 

Nos revenidos acima de 500ºC a martensite 
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priedades mecânicas do grupo de aços com 
12-15"/o Cr serão analisadas em 2.4. 

Se se desejar utilizar o aço na condição 
martensítica, um ciclo de tratamentos térmicos 
aconselhável é o seguinte: 


— Têmpera de 1050-1100 ºC, que permite obter 
martensite com todo o carbono em solução 
e, ao mesmo tempo, um tamanho de grão 
relativamente grande. 

— Têmpera de 900-950ºC, cuja finalidade é 
obter um refinamento de grão, mantendo 
a estrutura integralmente martensítica. 
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Fig. 3 — Variação da resistência à tracção (A), tensão a 0,1 "/; de deformação plástica 

(B), taxa de encruamento (C) e alongamento (D), com o tempo de tratamento a 500 “€ 

de um aço com 0,2 º, Ce 11,5º/ Cr. As faixas indicam os tempos ao fim dos quais pre- 
cipitam os carbonetos Cr;C; e Cr;;C, (dos resultados de [ 5]). 


transforma-se total ou parcialmente em sorbite, 
de acordo com a temperatura. À zona de amacia- 
mento máximo é bastante larga — 600 a 750ºC, 
sensivelmente — o que é uma vantagem sobre os 
aços de carbono, no aspecto de elaboração indus- 
trial, Neste intervalo de temperaturas o efeito 
maior é o crescimento das partículas de carboneto, 
com o aumento da temperatura e'ou do tempo 
de estacionamento. 

A escolha das temperaturas de tratamento, 
depende, evidentemente, das características estru- 
turais e mecânicas que se deseje obter. As pro- 
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— Revenido a 250ºC, para eliminação de tensões 
residuais de têmpera, 


A têmpera de alta temperatura possibilita a 
dissolução mais rápida do carbono e a homoge- 
neização da estrutura. A têmpera de 900-950 ºC 
elimina o inconveniente do tamanho de grão 
excessivo. O revenido a baixa temperatura poderá 
dispensar-se em alguns casos. Não será necessário 
naquelas peças que, pelo seu formato simples e 
pequena espessura, não estejam sujeitas a grandes 
gradientes térmicos. 
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Com este ciclo de tratamentos térmicos as ca- 
racterísticas mecânicas são as melhores, excepto 
em relação ao alongamento (ver 2.4). 


2.2 — Aços com teor em carbono inferior 
a 0,35 0h 


Com teores em carbono inferiores a 0,3-0,35 /o 
a estrutura de equilíbrio à temperatura ambiente 
é constituída por ferrite e perlite. Para teores 
inferiores a 0,1"/9 C a perlite tende a desapa- 
recer na sua forma característica — os carbone- 
tos tendem a formar grânulos, que aparecem 
rodeando os grãos de ferrite. 

Por arrefecimento ao ar, desde a zona auste- 
nítica, estes aços ou são integralmente martensí- 
ticos ou contêm, para os teores muito baixos em 
carbono, alguma ferrite. De um modo geral rea- 
gem aos tratamentos térmicos de maneira idên- 
tica à dos aços com 0,3-0,35 “A €. 


2.3 — Aços com teor em carbono superior 
a 0,35 %% 


O aumento do teor em carbono acima de 
0,35 "/o afecta considerâvelmente as temperatu- 
ras de transformação. Com 12-15 “h Cre1%4C 
a zona austenítica é já bastante estreita e situa-se 
a cerca de 1200 “C, iniciando-se a fusão a 1250ºC., 

As propriedades inoxidáveis destes aços não 
são muito boas. A elevada quantidade de carbo- 
netos de crómio que se forma reduz o teor 
efectivo de crómio (crómio em solução) de modo 
apreciável e daí resulta perda de resistência 
à corrosão. 

Altos teores em carbono permitem estabilizar 
a austenite à temperatura ambiente. A fig. 4 
permite verificar que, com 13 */y Cr, começa 
a reter-se austenite com cerca de 0,5 “/, C e que 
acima de 0,7 “/o C esta é totalmente retida. 

O efeito da temperatura de têmpera e de reve- 
nidos posteriores sobre um aço com 1,01 */% € 
e 11, 8 “/o Cr. está indicado na fig. 5. A dureza 
após têmpera auménta com a temperatura de 
têmpera até 1000-1075ºC (curvas A e B). Entre 
1100 a 1200 C (curva C) a dureza baixa, por 
retenção da austenite. Sobre as amostras marten- 
síticas — À e B— os revenidos têm efeito idêntico 
ao assinalado para conteúdos menores em car- 
bono. A amostra C, austenítica à temperatura 
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ambiente, mantém-se como tal até sensivelmente 
500ºC. Entre 500 e 580 parte da austenite 
transforma-se em martensite, que se mantém 
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Fig. 4 — Composição dos aços de crómio que forne- 
cem austenite à temperatura ambiente, quando arrefe- 
cidos rapidamente de 1200 *C [3]. 


à temperatura ambiente e provoca o aumento de 
dureza. O revenido a 600ºC conduz a troos- 
tite + martensite, e a 700ºC apenas a sorbite, 
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Fig. 5 — Efeito do revenido em um aço inoxidável de 
alto carbono, Arrefecimento ao ar de: A-g00 ºC; 
B- 1000 Ce C-1200 'C [3]. 


2.4 — Propriedades mecânicas dos aços com 
12-15 9/ Cr 


A tabela I indica os efeitos de diversas tem- 
peraturas de têmpera e revenido, para aços com 
composições comparáveis em crómio e carbono. 
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TABELA I 


Influência das temperaturas de têmpera e revenido nas características mecânicas 
de aços com 12-13 º/0 Cr [3] 


COMPOSIÇÃO 
(%) 


TRATAMENTO 
TÉRMICO 
TR á 
PEÇÃO exam?) 


RESISTÊNCIA A 


| 
0,22 060 013 0,14 


12.5 0,25 0.21 0,18 0,22 
| 


130 /016/0,31 0,52 /0,27 
| 


Para cada temperatura de revenido, o aumento 
da temperatura de têmpera conduz a maior resis- 
tência à tracção, menor alongamento e menor 
resiliência (comparem-se os casos a e b e os ca- 
sos c, d e e). O tamanho de grão y aumenta com 
a temperatura de têmpera e o material adquire 
maior temperabilidade. A estrutura resultante 
é, por isso, mais dura e menos dúctil. O aumento 
da temperatura de revenido, para a mesma tem- 
peratura de têmpera, produz efeito contrário: 
diminui a resistência à tracção e aumentam 
o alongamento e a resiliência (compare-se b 
com f). o 

As figuras 6 e 7 indicam as propriedades me- 
cânicas de aços com, respectivamente, 12-14 “/ 
Cr, C< 0,15 e 13-15 Cr, 0,26% C. 
A comparação das duas figuras permite verificar 
de imediato o efeito do aumento do teor em car- 
bono: a resistência à tracção e o limite elástico 
são mais altos, e o alongamento e resiliência 
mais baixos. Em ambas as figuras a curva de 
variações da resiliência com a temperatura de re- 
venido apresenta um pico característico a 250ºC, 
que não corresponde a qualquer flutuação nas 
curvas das outras características mecânicas. Este 
pico resulta da eliminação das tensões residuais 
de têmpera. Para revenidos a temperaturas supe- 
riores não deveria existir decréscimo de resiliên- 
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Fig. 6 — Características mecânicas de aços inoxidáveis 
martensíticos e sua variação com a temperatura de reve- 
nido (após têmpera em óleo de 1000 ºC) [6]. 
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